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I present a complete set of gauge invariant observables, in the context of general relativity coupled 
with a minimal amount of realistic matter (four particles) . These observables have a straightforward 
and realistic physical interpretation. In fact, the technology to measure them is realized by the 
Global Positioning System: they are defined by the physical reference system determined by GPS 
readings. The components of the metric tensor in this physical reference system are gauge invariant 
quantities and, remarkably, their evolution equations are local. 
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In general relativity (GR), the "observables", namely 
the physical quantities that we can predict and measure 
in real experiments, correspond to quantities of the the- 
ory that are invariant under coordinate transformations. 
It is not difhcult to construct observables in concrete ap- 
plications of the theory. Indeed, each physically mean- 
ingful number that we actually measure corresponds to 
one of these observables. For instance, in studying the 
general relativistic dynamics of the solar system, the dis- 
tance of a planet from us, at given (solar) time (say, at 
midnight tonight) is an observable. On the other hand, 
most of the theoretical work is not done in terms of ob- 
servables, but rather in terms of gauge dependent, that 
is, coordinate dependent, quantities. This of course is 
the reason for which the same physical situation can be 
described in terms of different metrics tensors. The usual 
procedure in GR is indeed to develop the theoretical de- 
scription of a certain physical situation in an arbitrary 
gauge, and then compute the value of a certain num- 
ber of coordinate independent observables, which can be 
compared with observations. 

This way of proceeding has an unsatisfactory aspect: 
we do not have at our disposal a complete set of observ- 
ables, which we could imagine, in principle, to measure 
in a simple manner. A physical situation, of course, is 
characterized by an equivalence class under diffeomor- 
phisms of solutions of the equations of motion; but we 
do not know how to effectively coordinatize the space of 
the equivalence classes in terms of realistically observable 
quantities. The individual observables that we use in con- 
crete applications are just a small number, and are far 
from capturing the full gauge invariant physics. The dif- 
ficulty of writing a complete class of observables is in fact 
well known (see [Q-^, and references therein), and raises 
several well known problems. For instance it complicates 
the canonical analysis of the theory and it is a serious 
obstacle to quantization. Furthermore, it generates con- 
ceptual difficulties for the very physical interpretation of 
the theory Q . 

Attempts to define a complete class of observables 
abound in the literature. It is easier to construct such 
observables in the presence of matter that in the context 



of pure GR. This is because the difficulty of writing ob- 
servables is the consequence of the absence of absolute 
localization in a general relativistic theory. If there is 
matter we can localize things with respect to the matter. 
For instance, we can consider GR interacting with four 
scalar matter fields. Assume that the configuration of 
the fields is sufficiently nondegenerate. Then the compo- 
nents of metric field at points defined by given values of 
the matter fields are observable. The idea is clearly the 
same as defining an observable as the distance between 
us and a planet: instead of having just a few planets to 
define a few physical distances, we imagine to have a con- 
tinuous of matter so that physical distances are defined 
everywhere. This idea has been developed in a num- 
ber of variants, such as dust carrying clocks and others, 
by many authors, including the present one. See 
and references therein. One succeeds in constructing a 
complete set of observables, but the extent to which the 
result is realistic or useful is certainly questionable. It is 
rather unsatisfactory to understand the theory in terms 
of fields that do not exist, or phcnomenological objects 
such as dust, and it is questionable whether these proce- 
dures could make sense in the quantum theory, where the 
aim is to describe Planck scale. Other (earlier) attempts 
to write a complete set of observables are in the con- 
text of pure GR |^] . The idea is to construct four scalar 
functions of the metric (say, scalar polynomials of the 
curvature), and use these to localize points. The value 
of a fifth scalar function in a point where the four scalar 
functions have a given value is an observable. This works, 
but the result is mathematically very intricate and phys- 
ically extremely unrealistic. It is certainly possible, in 
principle, to construct detectors of such observables, but 
I doubt any experimenter would get founded for propos- 
ing to build such apparata. Do we thus have to declare 
defeat and, in a general relativistic context, give up the 
hope of having a complete set of well defined observables, 
which are realistic, easy to construct, and do not assume 
that the world is different from what it is? 

In this paper, I'd like to argue that we do not have 
to declare defeat. I propose a simple way out, based on 
GR coupled with a minimal and very realistic amount 
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of additional matter. Indeed, the solution I discuss is so 
realistic that it is in fact real: it was inspired by an al- 
ready existing technology, the Global Positioning System 
(GPS), the first technological application of GR, or the 
first large scale technology that needs to take GR effects 
into account . The other source of the ideas presented 
here is the Null Surface Formulation (NSF) approach to 
GR §. What follows can be seen as a sort of pedestrian 
version of some NSF ideas. Ideas related to the ones pre- 
sented here have been presented by Massimo Pauri and 
Luca Lusanna ||] . See also . 

The idea here is simple. Consider a general covari- 
ant system formed eucaby GR coupled with four small 
bodies. These are taken to have negligible mass; they 
will be considered as point particles for simplicity, and 
called "satellites" . Assume that the four satellites follow 
timelike geodesies; that these geodesies meet in a com- 
mon (starting) point O; and, that at O they have a given 
(fixed) speed -the same for the four- and directions as the 
four vertices of a tetrahedron. The theory might include 
any other matter. Then (there is a region TZ of space- 
time for which) we can uniquely associate four numbers 
s",a = 1,2,3,4 to each spacetime point p as follows. 
Consider the past lightcone of p. This will (generically) 
intersect the four geodesies in four points pa ■ The num- 
bers s" are defined as the distance between pa and O. 
We can use the s"'s as physically defined coordinates. 
The components 5^/3 (s) of the metric tensor in these 
coordinates are observable quantities. They are invari- 
ant under four-dimensional diffeomorphisms (because, of 
course, these deform the metric as well as the satellites' 
worldlines). They define a complete set of observables 
for the region TZ. 

The physical picture is simple, and its realism is trans- 
parent. Imagine that the four "satellites" are in fact 
satellites, each carrying a clock that measures the proper 
time along its trajectory, starting at the meeting point O. 
Imagine also that each satellite broadcasts its local time 
with a radio signal. Suppose I am at the point p and have 
an electronic device that simply receives the four signals 
and displays the four readings. See Figure 1. These four 
numbers are precisely the four physical coordinates s" 
defined above. Current technology permits to perform 
these measurements with accuracy well within the rela- 
tivistic regime . If I then use a rod, and a clock and 
measure the physical distances between s" coordinates 
points, I am directly measuring the components of the 
metric tensor in the physical coordinate system. In the 
terminology of Ref. the s"'s are partial observables, 
while 30/3(5) are complete observables. 

As shown below, the physical coordinates s" have nice 
geometrical properties; they are characterized by 



g""(s) = 0, 



1, 



(1) 



Surprisingly, in spite of the fact that they are defined 
by what looks like a rather nonlocal procedure, the evo- 
lution equations for Qapis) are local. These evolution 



equations can be written explicitly using the Arnowitt- 
Deser-Misner (ADM) variables |l^] . Lapse and Shift turn 
out to be fixed local functions of the three metric. 




FIG. 1. si and S2 are the GPS coordinates of the point p. E 
is a Cauchy surface with p in its future domain of dependence. 

In what follow, we begin for simplicity by introducing 
the GPS coordinates in Minkowski space. This allows 
us to introduce some tools in a simple context. Then 
we go over to a general spacetime. We take the speed 
of light to be one, signature — , — , — ], and we assume 
the Einstein summation convention only for couples of 
repeated indices that are one up and one down. Thus a 
is not summed over in (|]) . While dealing with Minkowski 
spacetime, the spacetime indices /i, v are raised and low- 
ered with the Minkowski metric. We write an arrow over 
three- as well as four-dimensional vectors. 

Consider a tetrahedron in three-dimensional euclidean 
space. Let its center be at the origin and it four ver- 
tices t/", where and v"- ■ = —l/i for a ^ (3, have unit 
length \v"-\'^ = 1. Here a = 1, 2, 3, 4 is an index that dis- 
tinguishes the four vertices, and should not be confused 
with vector indices. With a convenient orientation, these 
vertices have polar coordinates u = (r, 9, (f) 



^^'" = (1,0,0) 
w^" = (l,a,0) 
w^" = (l,a,27r/3) 



V 



(l,a,-27r/3) 



with cos a 
1,2,3) 



-1/3, and cartesian coordinates 



,1a 



(0, 0, 1) 
v'^" = (2^/2/3, 0, -1/3) 

= (-V2/3, v/273, -1/3) 
v^"- = (-\/2/3, -v/273, -1/3). 
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Let us now go to a four dimensional Minkowski space. 
Consider four timelike 4- vectors W", of length one, 
|yyQ!|2 _ representing the normalized 4- velocities of 
four particles moving away from the origin in the direc- 
tions i/" at a common speed v. Their Minkowski coordi- 
nates (/X — 0, 1, 2, 3) are 
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(1, 



(10) 



We fix the velocity v by requiring the determinant of the 
matrix W"^ to be one one. (This choice fixes v at about 
half the speed of light; a different choice changes only a 
few normalization factors in what follows.) The four by 
four matrix W"''^ plays an important role in what follows. 
Notice that it is a fixed matrix whose entries are certain 
given numbers. 

Consider one of the four 4- vectors, say W = . Con- 
sider a free particle in Minkowski space that starts from 
the origin with 4- velocity W . Call it a "satellite". Its 
world line I is a;(s) = sW . Since W is normalized, s is 
precisely the proper time along the world line. Consider 
now an arbitrary point p in Minkowski spacetime, with 
coordinates X. We want to compute the value of s at 
the intersection between I and the past light cone of p. 
The calculation is particularly easy in a Lorentz frame 
in which the satellite stays still, namely has vanishing 
3-velocity, and p is in the (t, x) plane, with coordinates 
(T, X). Assume X > 0. In this frame, the equation of I 



a; = 0, 



(11) 



the proper time is s = t, and the equation of the (inter- 
section with the (t, x) plane of the) past light cone of p 
is 

(X-a;) = ±(r-t), t<T. (12) 
Taking all this together gives 

s = T-X. (13) 

Notice now that in this frame |Xp = T^-X^ and X-W = 
T, so that we can write 



T - X = T - vr23(r23^ 



= X-W-^{X-Wf-\X\^. (14) 
From the last two equations we obtain 



s = X - W - d {X - Wf ~\X\ 



(15) 



But this equation is Lorentz covariant and therefore it is 
true in any Lorentz system. (A direct full four dimen- 
sional calculation in an arbitrary system gives the same 
result.) 



Let us now consider four satellites, moving out of the 
origin at 4- velocity W^". If they radio broadcast their 
position, an observer at the point p with Minkowski co- 
ordinates X receives the four signals s" 



^X-W 



{x-w°'Y-\x\ 



(16) 



We introduce (non-Lorentzian) general coordinates s" on 
Minkowski space, defined by the change of variables (|l^). 
These are the coordinates read out by a GPS device in 
Minkowski space. The Jacobian matrix of the change of 
coordinates is given by 



dxf^ " 



W^{X -W") - X^ 



{x-w)^~\x\-^ 



(17) 



where and X^ are W"^ and X^ with the spacetime 
index lowered with the Minkowski metric. This defines 
the the (co-)tetrad field £^^(s) 



(18) 



The contravariant metric tensor is given by g"^{s) = 
E^{s)Ei'f^{s). Using \W°'\'^ = 1, a straightforward calcu- 



lation shows that 



5""(s) = 0, a-l,...,4. 



(19) 



This equation has the following nice geometrical interpre- 
tation. Fix a and consider the one- form field = ds". 
In s" coordinates, this one-form has components oj'^ — 
S^, and therefore "length" jcj^p = g^Tw^w^ = g°"^. 
But the "length" of a 1-form is proportional to the 
volume of the (infinitesimal, now) 3-surface defined by 
the form. The 3-surface defined by ds" is the surface 
s" — constant. But s" = constant is the set of points 
that read the GPS coordinate s", namely that receive 
a radio broadcasting from a same event pa of the satel- 
lite a, namely that are on the future light cone of pa- 
Therefore s" — constant is a portion of this light cone, 
therefore it is a null surface, therefore its volume is zero, 
therefore = 0, therefore 5"" = 0. 

Since the s" coordinates define s" = constant surfaces 
that are null, we denote them as "null GPS coordinates" . 
It is useful to introduce another set of GPS coordinates 
as well, which have the traditional timelike and spacelike 
character. We denote these as s**, call them "timelike 
GPS coordinates" , and define them by 



(20) 



This is a simple algebraic relabeling of the names of the 
four GPS coordinates, such that s'^=° is timelike and 
gtj.=a spacelike. In these coordinates, the gauge condi- 
tion (|l9| ) reads 



W^W^g^^''{s)^0. 



(21) 
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This can be interpreted geometrically as follows. The 
(timelike) GPS coordinates are coordinates such that 
the four 1-forms fields 



(22) 



are null. 

Let us now jump from Minkowski space to full GR. 
Consider GR coupled with four satellites of negligi- 
ble mass that move geodesically and whose world lines 
emerge from a point O with directions and velocity as 
above. Locally around O the metric can be taken to be 
Minkowskian; therefore the details of the initial condi- 
tions of the satellites worldlines can be taken as above. 
The phase space of this system is the one of pure GR plus 
10 parameters, giving the location of O and the Lorentz 
orientation of the initial tetrahedron of velocities. The 
integration of the satellites' geodesies and of the light 
cones can be arbitrarily complicated in an arbitrary met- 
ric. However, if the metric is sufficiently regular, there 
will still be a region TZ in which the radio signals broad- 
casted by the satellites are received. (In case of multiple 
reception, the strongest one can be selected. That is, if 
the past light cone of p intersects I more than once, gener- 
ically there will be one intersection which is at shorter lu- 
minosity distance.) Thus, we still have well defined phys- 
ical coordinates s" on TZ. Equation ( pj| ) holds in these 
coordinates, because it depends just on the properties of 
the light propagation around p. We define also timelike 
GPS coordinates by (20), and we get the condition 
( pT| ) on the metric tensor. 

To study the evolution of the metric tensor in GPS co- 
ordinates it is easier to shift to ADM variables N, N"-, jab- 
These are functions of the covariant components of the 
metric tensor, defined in general by 

ds^ — g ^i,dx^ dx^ 

= N^dt^ - jabidx" - N''dt){dx^ - N^dt). (23) 

Equivalently, they are related to the contravariant com- 
ponents of the metric tensor by 



(24) 



where 7°'' is the inverse of jab and — (1/iV, N°-/N). 
Using these variables, the gauge condition (eTJ) reads 



(W"n^)2. 



(25) 



Notice now that his can be solved for the Lapse and Shift 
as a function of the three-metric (recall that are fixed 
numbers), obtaining 



Or, exphcitly. 



N : 



(28) 



The geometrical interpretation is as follows. We want 
the 1-form cj" defined in (^2|) to be null, namely its norm 
to vanish. But in the ADM formalism this norm is the 
sum of two parts: the norm of the pull back of on 
the constant time ADM surface, which is q", given in 
(^), and depends on the three metric; plus the square 
of the projection of uj" on n^. We can thus obtain the 
vanishing of the norm by adjusting the Lapse and Shift. 
We have four conditions (one per each a) and we can 
thus determine Lapse and Shift out of three metric. In 
other words, whatever is the three metric, we can always 
adjust Lapse and Shift so that the gauge condition ( [2l| ) 
is satisfied. 

But in the ADM formalism, the arbitrariness of the 
evolution in the Einstein equations is entirely captured 
by the freedom in choosing Lapse and Shift. Since here 
Lapse and Shift are uniquely determined by the three 
metric, evolution is determined uniquely if the initial 
data on a Cauchy surface are known. Therefore the evo- 
lution in the GPS coordinate s° of the GPS components 
of the metric tensor, g^v{s), is governed by deterministic 
equations: the ADM evolution equation with Lapse and 
Shift determined by equations (p7|-p8|) . Notice also that 
evolution is local, since the ADM evolution equations, as 
well as the equations (p7|-|2^), are local. [| 

To understand what is going on in this gauge, consider 
an arbitrary coordinate system x^ , and an arbitrary co- 
ordinate transformation s" = 3 " (x) such that the trans- 
formed metric tensor satisfies (19). This is equivalent 
to imposing one differential equation on each of the four 
functions s"(a:), as follows 



g^"'{x)d^,s°'{x)d,,s°'{x) = 0. 



(29) 



Equivalently, the problem is to find four independent ex- 
act null 1-forms = ds", integrable over a finite re- 
gion. Or, equivalently again, the problem is to find four 
nowhere-parallel foliations of Minkowski space with fam- 
ilies of null surfaces. The last is precisely the main in- 
gredient of the Null Surface Formulation of GR 0. A 10 
parameter family of solutions of this problem is given by 
the construction above, namely by the light fronts emerg- 
ing from the geodesies of the satellites, and labelled with 
the proper time along these geodesies. Equation (p^ ) has 
of course more solutions than these, since there are more 
families of null foliations than the ones constructed here. 
A null surface in a solution of (E9), for instance, does 



Wiiq'^ 



(26) 



where Wj^ is the inverse of the matrix and 



(27) 



*This does not imply that the full set of equations satisfied 
by (;^ii/(s) must local, since initial conditions on s° = satisfy 
four other constraints besides the ADM ones. 
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not need to focus in a point. Even less, the surfaces of 
the family need to focus on a line which is a geodesic, 
and so on. But if we require that the surfaces focus on a 
point, that the points are along a timelike geodesic, that 
the label of the foliation is given by the proper time along 
this geodesic, and that the four geodesies meet at a single 
point and with the prescribed angles, then the space of 
solutions of ( p9|) is reduces down to a finite 10 parameter 
space. The 10 parameters are the location of the origin O 
and the Lorentz orientation of the tetrahedron of initial 
velocities. (By pushing O all the way to past infinity, one 
obtains the asymptotic NSF coordinates.) 

How can the evolution of the quantities giiv{s) be lo- 
cal? The conditions on the null surfaces described in 
the previous paragraph arc nonlocal. Coordinate dis- 
tances yield typically to nonlocality: Imagine we define 
physical coordinates in the solar system using the cos- 
mological time tc and the spatial distances xs,xe,x,j (at 
fixed ic) from, say, the Sun the Earth and Jupiter. The 
metric tensor g^v{tc,xs,XE,xj) in these coordinates is 
observable, but its evolution is highly non local. To see 
this, imagine that in this moment (in cosmological time), 
Jupiter is swept away by a huge comet. Then the value 
of 5/j!/(ic, xs,xe,xj) here changes instantaneously, with- 
out any local cause: the value of the coordinate xj has 
changed because of an event happened far away. What's 
special about the GPS coordinates that avoids this nonlo- 
cality? The answer is that the value of a GPS coordinate 
in a point p does in fact depend on what happens "far 
away" as well. Indeed, it depends on what happens to the 
satellite. However, it only depends on what happened to 
the satellite when it was broadcasting the signal received 
in p, and this is in the past of p ! If p is in the past do- 
main of dependence of a partial Cauchy surface S, then 
the value of g^v{s) in p is completely determined by the 
metric an its derivative on E, namely evolution is causal, 
because the entire information needed to set up the GPS 
coordinates is in the data in S. See Figure 1. Explicitly, 
the = constant surfaces around E can be uniquely 
integrated ahead all the way to p. They certainly can, as 
they represent just the evolution of a light front! This is 
how local evolution is achieved by these coordinates. 

Let us summarize. We have introduced a set of phys- 
ical coordinates, determined by certain material bodies. 
Geometrical quantities such as the components of the 
metric tensor expressed in physical coordinates are of 
course observable. This procedure is well known. The 
novelty here is that we have shown that in order to ob- 
tain physical coordinates, there is no need, as usually 
assumed, to introduce a large unrealistic amount of mat- 
ter or to construct complicated and unrealistic physical 
quantities out of the metric tensor. Indeed, four par- 
ticles are sufficient to coordinatize a (region of a) four- 
geometry. Furthermore, the coordinatization procedure 
is not artificial: it is the real one utilized by existing 
technology. 



Notice that the degrees of freedom of the theory we 
have considered are still only two per point (plus the 10 
parameters of the initial condition of the satellites.) This 
is why initial data on for g^v{s) must satisfy additional 
constraints. We leave the study of these and the explicit 
construction of the associated Hamiltonian formalism for 
future investigations. 

We conclude with some simple comments on observ- 
ability. The components of the metric tensor in (time- 
like) GPS coordinates can be measured in principle as 
follows. Take a rod of physical length L (small with re- 
spect to the distance along which the gravitational field 
changes significatively) with two GPS devices at its ends 
(reading timelike GPS coordinates). Orient the rod (or 
search among recorded readings) so that the two GPS 
devices have the same reading s of all coordinates except 
for s^. Let 5s^ be the difference in the two readings. 
Then we have along the rod 

ds^ = gii{s)5s^5s^ ^ L^. (30) 

Therefore 

g^^{s) = {L/5s'f. (31) 

Non-diagonal components of gab{s) can be measured by 
simple generalizations of this procedure. The gobis) are 
then algebraically determined by the gauge conditions. 
In a thought experiment, a spaceship could travel in a 
spacetime region and compose a map of values of the 
GPS components of the metric tensor. 

To avoid confusion, it is may be useful to recall here 
the distinction between partial and complete observables 
]lT[ . A partial observable is a quantity to which a mea- 
suring procedure can be associated. A complete observ- 
able is an observable quantity that can predicted by the 
theory, or, equivalently, whose knowledge provides infor- 
mation on the state of the system. For instance, in the 
theory of a single harmonic oscillator q{t) the time t and 
the position of the oscillator q are partial observables, 
while the position of the oscillator at a given time, or at 
the initial time, is a complete observable. The GPS coor- 
dinates are partial observables: we can associate them a 
measuring procedure (this is what has been done in this 
paper), but we can of course not "predict" s^. Equiv- 
alently, there is no sense in which the metric field of 
the universe is characterized by a certain value of s^. 
The complete observables, or true observables, are the 
quantities g^ivis), for any given value of the coordinates 
s^. These quantities are diffeomorphism invariant, are 
uniquely determined by the initial data and in a canoni- 
cal formulation are represented by functions on the phase 
space that commute with all constraints. 

Finally, notice that the observables we have defined 
are a straightforward generalization of Einstein's "point 
coincidences"^]. In a sense, they are precisely Einstein's 



"All our space-time verifications invariably amount to a 
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point coincidences. Einstein's "material points" are just 
replaced by photons (light pulses): the spacetime point 
s" is characterized as the meeting point of four photons 
designated by the fact of carrying the radio signals s". 
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